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THE ROLE OF ATOMIC OXYGEN IN THE IONOSPHERIC 
E AND F REGION BEHAVIOR 

J. R. Herman and S. Chandra 
Laboratory for %>ace Sciences 

ABSTRACT 

The continuity and heat conduction equations for the elytron, ion* and 
neutral gases are solved simultaneously to show the effect of using the heat loss 

e 

arising from the excitatim of the fine structure levels of atomic m^gen, and the 
heat iiqmt from a height dependent electrcm heating efficiency. The solutions 
are compared to Thomson backscatter data obtained under conditi(»s similar to 
those used in the theoretical model. A family of soluticais is presented for three 
di!ierent values of the atomic oi^gen density at the turbopause level, and com- 
pares ffivorably to certain widely observed features of the main phase of a mag- 
n^ic storm. From the comparison, it is shown that the results of decreasing 
the atomic oi^gen density at the turbopause level are consistent with the ob- 
served changes in both the neutral and iooized constituents relative to quiet 
magnetiw conditions. 
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THE ROLE OF ATOMIC OXYGEN ON THE IONOSPHERIC 


E AND F REGION BEHAVIOR 

INTRODUCTION 

In a recent paper, Dalgarno and Degges (1968) have shown that an efficient 
way of coolii^ the electron gas in the E and F regions of the ionosphere is by 
the excitation of the firo structure levels of atomic oxygen through inelastic col- 
lisions with thermal electrons. This loss process is ccmsiderably more impor- 
tant than other ekistic and inelastic mechanisms usually considered for con- 
struction of theoretical models of the electron temperature. The use of this ad- 
ditional loss term in the electron heat conduction equation offers the possibility 
of accounting for some of the discrepancies which exist between the theoretically 
derived electron temperature profiles and those measvired by techniques such as 

Thomson backs catter (Dalgamo et al., 1967 and 1968). 

Recently, in presenting the results of the simultaneous solution of the elec- 
tron continuity equation and the heat conduction equations for the electron, ion, and 
neutral gases, we briefly discussed the effect of cooling the electron gas thmngh 
excitation of the atomic oxygen fine structure levels (to be abbreviated here- 
after as A.O.F.S. loss). Our result showed that this dominant cooling mechanism 
has the effect of considerably reducing the electron temperature in the 150 to 
250 km altitude range and to a lesser extent elsewhere (Herman and Chandra, 
1969; to be referred to as paper Q. The purpose of this work is to show the ef- 
fects of including the A.O.F.S. loss process in the coupled equations for T^, T. , 
and ri^, and to discuss some of the implications for the general behavior of the 
F-region. 
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The Equations of Continuity and Heat Conduction 

The continuity and heat conduction equations, their boundary conditions, and 
the method of obtaining steady state solutions as applied to the F-region ionos- 
phere have been discussed in detail in paper 1. For convenience an outline of 
these equations is presented using the same notation and numerical values. 



These equations are the electron continuity equation, the heat conduction equa- 
tions for electron, ion, and neutral gas temperatures, and the equation of diffusive 
equilibrium for the neutral gas densities, respectively. 


2 


The new term in Equation 2, L^* , represents the A..O.F.S. loss, and is given 
by the expression (see appendix) 
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= die electron density 
= the atomic oxygen density 
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, T. , are the electron, ; and neutral temperatures 
k = Boltzmann's constant 

Ej = the i*** fine structure energy level of atomic oxygen cor- 
re^onding to the angular momentum quantum number j 
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gj = 2j + 1 

^eo ( i ' ** i > E' ) = the cross section for downward transitions from j' to j 

D=2^e,.e:v(-E|,/kT„) 

j' 


Equation 6 has been obtained by combining the terms for upward and downward 
transitions using detailed balance and then averaging over the Maxwellian velocity 
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distribution of the electrons. In addition, it is assumed that the population of 

atomic oxygen fine structure energy levels is given by the Boltzmann distribu- 

« 

tion at temperature T^. When the numerical values for ( j ' - j . E' ), given 
by Breig and Linn (1966), are used, the following numerical approximation to (6) 
is obtained (see appendix) 

^ (S. 56 - 2. 40 X 10'* T ) ((} . 23 ^ 2. 75 x 10 • , ) 

X 10*^®n^N^(T^ )/^n ^*"85 cm‘^ sec”^ T(7) 

which agrees approximately with the values given *>y Dalgarno and Degges (1968). 

For Equations 1 to 5 the bnating ii cd ion elation efficiencies are available as 
free parameters if adjustments are t.;;eded lor comparison with data. 'Vith the 
use of Equation 7 it was felt that the ionization efficiency should be taken as 1.0 
(see paper I) along with two cases of the electron heating efficiency. The heat 
production function, Q^, is obtained by computing the production rate of eiectron- 
ion pairs and multiplying by a constant (with altitude) heating efficiency of 2.5 ev 
per electron ion pair. Since this procedure underestimates the heating rate in 
the high altitude range, a height dependent, electron heating efi^iciency suggested 
by Dalgarno et al. (1968), and approximated by the values shown in Figure 1, is 
used for comparison with the constmit efficiency. 

DISCUSSION 

The results of solving Equations 1 to 5 with and without the A.O.F.S. loss 
and the height dependent electrcm heating efficiency are shown in Figure 2. The 
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parameters used for this case correspond to low solar activity (HF* = G.9);^ 
afternoon conditions (time = 13:30 hours), latitude and dip approximating 
Arecibo, Puerto Rico (Lat = 18® and Dip = 50®), and for the solar declination 
angle computed for September 22. In addition, the boundary conditions for the 
neutral gas density at 100 km were taken as [Nj] = 4.65 10^^ , [O^j = 1,25 x 10*^ 

[O] = 7.50 xio“, [^e] = 2.50 X 10®, and [H] = 1.00 x lo® cm"®. 

The most significant change occurs in the vicinity of 190 km where the high 
peak electron temperature is reduced by 50%. When only the A.O.F.S. loss is 
used, the significant temperature changes are confined to the lower altitudes 
below about 300 km as is shown in Figure 2 (curve 2) . However, with the height 
dependent electron heating efficiency and the A.O.F.S. loss, the high altitude 
electron temperature is enhanced by about 30% (see curve 3). The temperature 
enhancement is in turn reflected in a considerable change in the slope of the 
electron density profile above F2 . The higher temperature in the region where 
diffusion is important produces the change in shown in the upper portion of 
curve 3. The large' decrease in below 300 km produces only a small change 
in n^, since here the electron density profile is largely determined by chemical 
processes and not diffusion. The introduction of the A.O.F.S, loss does shift 
F2 downward by about 30 km with a slight increase in F2 , and also causes 
a decrease in layer thickness. 

‘The abbreviation HF* " N stands for the solar flux values in the spectral range of 10 to 1027 A, 
as given by Hinterregcr et al. (1965), tnultiplied by N. 



Since the heat input from the electron to the neutral gas amounts to about 
15% of the total, it mi^t be expected tiiat the addition of the A.O.F.S. loss would 
add to this heat input. The neutral gas temperatures corresponding to the curves 
of Figure 2 do not show such an effect (t^ (1) = 887®K, = 870*K 

since the heat input from the electron gas below about 300 km 
is largely determined by the magnitude of Q^. That is, with the addition of the 
A.O.F.S. loss the electron temperature adjusts itself downward until the heat 
loss rate to the other gases is about the same as before (Q^ is almost uneffected 
by A.b.F.Sv loss). Above about 400 km, where is not controlled by the A.O.F.S. 
loss, the electron ^tnperatiire- shews only a very small increase due mainly to 
the small decrease in electrvu density and neutral temperature (compare curves 
1 and 2). Th’S is horne out curve 3 (t^ (3) - 925°k), where the increased 
heat input at almost all altitudes, compared to curve 2, results in a high'^r 
value of Tg and in spite of increased n^. , 

For quantitative comparison of the present theoretical models with observa- 
tional data, it is necessary to simultaneously measure the density and tempera- 
ture of both ionized and neutral constituents along w ith the solar EUV flux. In 
the case of Thomson backscatter, the quantities that are usually measured are 
the charged particle temperatures and densities. The information at ut the 
solar EUV flux and the neutral gas composition have to be inferred independently 
from solar decimeter flux, satellite drag, and atmospheric models. With this 
limitation in view, it is meaningful to compare the theoretical results with the 
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Thomson backscatter data whenever the conditions during the eiqmriment are 
reasonably close to those used in the present calculation. An example of equinoc- 
tial midday data taken during a time of low solar activity ^pt. 22, 1965 at 13:15) 
has been presented fay Mahajan* (1967), and is reproduced in Figure 3. A com- 
parison witti the curves labelled 3 in Figure 2 (redrawn in Figure 3) shows a suf- 
ficient similarity to indicate that the introductimi of the A.O.F JS. loss and the 
height dependent electron heating efficieimy are nacessary for olAaining satis- 
hictory iono^dieric models. 

Another example simultaneous electitm density and temperature meas- 
urements from Thomson backsf^atter has been given by Rao (1968) for July 28, 
1965 at 15:00 hours (not reproduced here). Ibe conditions for this data rou^ly 
corre^^nd to the same solar zenith angle and solar activity as adopted in Fig- 
ure 2. hi this and other similar data taken under conditions that are well enou^ 
known to be compared to tte models, Uie main features are well reproduced for 
die T^, Tj, and n^ profiles. In addition to the magnitudes of T^, Tj. and n^, both 
the data and the theoretical models show a decrease in the logarithmic slope 
with increasing altitude whenever there are substantial electron temperature 
gradients above b„ F2. A similar change in the logarithmic slope can oo^ur 
when there are changes in the ionic composition from O* to H However, it is 
reasonable to assume that the ionic composition was mainly until above about 
750 km for the data taken on July 28, 1965 >Iuring the afternoon, and therefore 

*The electron density profile shown in Figure 3 is unpublished, and has been kindly prov ided by 
the author for comparison. 
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the decrease in the magnitude of the logarithmic slope can be attributed to the 
electron temperature gradient. " ~ 

THE EFFECTS OF VARIABLE [O] 

The effects of changing the atomic o^^gen density boundary values have been 
discussed previously in the absence of A.C F.S. loss and the hei^t d^endent 
electron heating efficiency (paper I also ClKindra and Herman, 1969; to be referred 
to as paper IQ. With the introduction of the A.Q.F.S. loss, both the neutral and 
electrcm gas temperatures aie directly controlled by the atomic oxygen density. 
Changes in [Ol, either throu^ dynamic or chemical processes, can be expected 
to have profound effects (m the E and F-region iono^here 

As ejqiected, the use of the A.O.F.S. loss and the hei^t dependent electron 
heating efficiency did not change any of the essential conclusions presented in 
either paper I or paper If, even though the details of the , T. « T^. and pro- 
files are substantially different. Figure 4 ^ows a family of profiles obtained 
by setting the atomic oxygen doisity at 100 km to 1,00 10**, 7.50 * 10** , and 

5.00 10** for curves A, B, and C rei^ctively. The remaining parameters 

have been held constant (for this case the pairameters correspond to low solar 
activity (HF* = 1.00), midday ccmditions (time = 12:00 hours), latitude and dip 
corre^ondii^ to Wallops Island, Virginia (Lat = ZT and Dip = 70®), and for the 
solar declination an^e computed for March 2. 

If, as discussed in paper n, the conditiontS well after the onset of a magn'^tic 
storm are chairacterized by a decrease in the atomic oi^'gen density at the 
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turbopause level, then curves A, B, and C can be interpreted as corresponding 
to quiet, disturbed, and very disturbed conditions, respectively. The commonly 
known changes in the F-regi<m parameters in the middle and high latitudes as- 
sociated with the main phase of a magnetic storm are given in the following list. 

1. A decrease in foF2 (tq> to 50%). 

2. An increase in F2 (up to 100 kin) . 

3. An increase in the topside electron density from about 100 km above the 
F2peak. 

4. A decrease in the total clectrcui content. 

5. A worldwide increase in satellite drag that can be interpreted as an in- 
crease in both neutral temperature and density. 

These features of a magnetic storm, as referred to quiet conditions (curve A), 
are clearly shown in Figure 4. In addition, the exospheric neutral temperature 
increased by over 200’K (t„ (A) = 796‘K, T„ (B) = 888*K, and T„ (C) = 1000*K) 
accompanied by an increase in the tc^al neutral density at-all heights. Also, as 
would be expected, the effect of decreasing atomic oxygen is to increase both the 
electron and i<m temperatures and the hei^t of the maximum. The main dif- 
ference between this and the previous descr4>tion of a magnetic storm by 
Chandra and Herman (1969) is in the magnitudes of the resulting temperatures 
and densities. %>ecifically, the rarely observed sharp peaks in the electron 
temperature no longer ai^ar, the gradients in both the electron temperature 
and electron density above 300 km increase, and the ion temperatures are 
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\ 

c(»siderably increased above 450 km where T. begins to tend stron^y 
towards and away from (note, however, that T. does not become equal 
to for the altitude range considered, 100 to 1000 km). 

SUMMARY 

The equations of heat conduction for the electron, ion, and neutral gases 
have been solved simultaneously with the electrcm continuity equati<m using the 
A.O.F.S. loss process and a height dependent heating efficiency. The resulting 
solutions show that the electron temperature is significantly reduced below about 
300 km where the A.O.FJS. loss is most important and increased at hif^er altitudes 
because of the increased heat input from the height dependent electron heating 
efficiency. Smaller changes take place in the electron density profile and the 
hei^t of the maximum, F2, because of the reduced electron temperature 
and gradients entering into the difrision velocity. Finally, the i;m temperature 
is increased in proportion to the increased electron temperature. 

The comparison of these solutions with Thoms(m backscatter data show that 
the changes brought on by including the A.O.FJ5. loss and the height dependent 
electron heating efficiency remove some of the qualitative and quantitative 
discrepancies between the theoretical and measured profiles. 

A family of profiles are presented for , T; , and for three different 
values of the neutral atomic oxygen density boundary condition at the turbopause 
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level. These results compare favorably with certain widely observed features 
L. the F- region ai^aring well after the onset of a magnetic storm. From 
this agreement it can be concluded that most of the observations can be eiqplained 
by a decrease in the atomic oxygen density at the turb<^use level r^ative to 
quiet magnetic conditims. 
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APPENDIX 

Consider a gas, each atom of which possesses internal energy eigenvalues 
interacting via. inelastic collisions with an electrou gas so as to cause transi- 
tions between the states E. and E. . Assume the gases to have kinetic energy 
distribution^ f , (e, ) and f , (Ea)» ci-'oss sections for the transitions between Ej 
and Ej of ( i - j , E) and ( i ** j . E) , for electron-atom and atom-atom col- 
lisions re^>ectively, and that the state E. has an occupation number n^. The 
number of transitions into the state i per unit volume per second is given by 


dn.. 

~3T 


- I d^ V id^ V n. f f 
I el a i e a 



(i-j) 



n. 

j 


f f 


e a 





Jd^ v.'JdS V^n. 


i ^ K) f. K')^. (i ^i) 



Jd3 V,' v, f , (v. ) f. (v; ) «i, ( j - i ) I V. - v; j 


(Al) 


Since cr( i -> j ) is assumed to depend only on the magnivude of the relative 
velocity, some of the integrals may be evaluated to yield 


dn.. 

"3T = 


-J dEf^ (E)n. (i-j.E) + J dE f ^ (E) |/^n. (j- i. E) 


f dEf, (E)Oj (i-j.E)+ [* dEf, (E)n. (j-i.E) 

Ji.j “ -0 ^ “ 


(A2) 
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where Ijj = E- -Ej, (E) is Maxwellian, | - v^' j - |^4E/mg, and and 

are the electron and atom masses respectively. 

During the inelastic collision only one of the particles is assumed to undergo 
an internal energy change white the other changes its kinetic energy from E to 
E' . Detailed balance can now be used to relate the upward and downward transi- 
tion cross sections. 

= E' g. <7(j - i, E') (A3) 

where g. is the degeneracy of the state. If it is assumed that E. > E., then 
the rate at which electrons lose energy to the atoms per unit volume is 

~^i) "3F (A4) 

i 

Ej>Ej 

Using A3 to combine the integrals in A2 by eliminating ( i j , E), assuming 
to be Maxwellian at a temperature and assumii^ that the electron colli- 
sicns are infrequent enou^ so that the occupation number n . can be well ap- 
proximated by the Boltzmann distribution 

i 8j «*xp ("Ej/kT. ) • (A5) 


14 


Equation A4 becomes 


l;. 


2 "eN, 






Ej>Ej 


1 - exp (t,, - T, 


. . ^ -E/kT 

Ea^ (j^i,E)e dE» 


(A6) 


(A6 is an alternate form to that given by Dalgamo and De^s, 1968), where 
and are the electron and atom munber densities respectively, and D is given 
by 


D = (-EjA'^'a) (A7) 

i 

When the numerical values for ( j - i , E) are used, as given by Breig and 
Linn (1966), the cross sections can be represented by the form 



E) 

= Aj, E» * E * C,, 


. (A8) 


E) 

-a 

o 


(A9) 

(i " 

E) 

= W„,ryE 

(EikTj 

(AlO) 


where kT^ represents an electron kinetic energy above which the cross sections 
are assumed to be a constant value or to decrease with a 1/E energy dependence. 
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i . 

! 

I The result of using Equations A8 and A9 to evaluate A6 is 

(®i ■ ®i) exp iTf;; - ft;) 

j i^j 
Bj>E. 

• ctV (T,-T.) 

\ e a 

F,j = [6A„(kT.)=*2B,j(kT.)*C„*(<7jx„*l)-A,,(kT.)'(x/‘3X/ 

t6X,t6)-B,j(kTj(x.2*2X„*2)-C,.(X,M)}*,’‘”] ,A12) 

K ■ W./tkTj 

If instead, AlO is used in place of A9, the term c.-^ replaces (X^ + l). For 
the case of A.O.F.S. loss, tiie values of A. ■ > B. . , and C. . are taken to be (for 
cr^ in units of the Bohr area, tts^, and E in ergs) 

= - 1.400x 102^ = 2.944 x10^2 q.3093 

Ajj = -9.171x1023 Bjj = 1.978 x10*2 = 0.5426 

Ajj, = - 8.814x 1023 Bj^ ^ 1.844x 10*2 = 0.2588 




Inserting these values in Equation All and performing the indicated summations 
for i = 1, 2, 3 and j = 1, 2, 3 give the following approximation for , 

Vo "" {5. 56 - 2. 40 X 10"* Tj (6. 23 + 2. 75 x 10" 3 T„) 

X 10'25 (T^ - T„ )/t„ ergs cm'^ sec"! (a13) 

When Equation AlO is used in place of A9, the same values (within 2'Je) of 1-/^ 
are obtained for electron temperatures less .tlmn 2500*’K. Above this temperature 
the results slowly diverge to produce values that are lower by as much as a 
factor of 1.8 when = 11,000"K. 
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Figure 1. 
Figure 2. 


Figure 3. 
Figure 4. 
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FIGURE CAPTIONS 

Height dependent electron heating efficiency. 

A comparison of constant heating efficiency without A.O.F.S. loss 
(curve 1) with constant heating efficiency with A.O.F.S. loss (curve 2) 
and height de^ndent electron heating efficiency with A.O.F.S. loss 
(curve 3). 

Thomson backscatter results for September 22, 1965 (13:15 hours) 
at Arecibo, Puerto Rico (Mahajan, 1967). 

The effect of varying the density of neutral atomic ojg^gen at the 
lower boundary (100 km). Curve A — [O] = 1.00 x 10^^, curve B — 
[O] = 7.50 X 10^^ , and curve C — [O] 5.00 x 1(7“ cm"^. 
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